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1.  INTRODUCTION 


This  is  one  of  a  series  of  reports  dealing  with  modeling  of  the  marine 
atmospheric  boundary  layer.  The  basics  of  the  tfeval  Rcstgraduate  School 
modeling  approach  are  given  in  Ehirall  (1981).  That  report  also  describes  the 
current  status  of  modeling  the  pertinent  physical  processes.  The  utility  of 
the  model  for  tactical  use  and  initial  model  validation  are  given  in  Davidson  et  al 
(1982).  This  report  is  in  a  samevAat  different  vein.  The  other  reports  deal 
with  use  of  the  model  to  describe  or  predict  the  evolution  of  boundary  la>er 
properties.  Hare  we  use  the  model  "in  reverse"  to  determine  an  important 
physical  parameter,  the  rate  of  generation  of  marine  aerosols . 

There  are  two  components  to  aerosols  in  the  marine  boundary  layer:  1)  con¬ 
tinental  (background)  and  2)  locally  generated  sea-spray  droplets  ( tomhardt 
and  Streete,  1970)  -  The  sea-spray  droplets  are  generated  primarily  by  the 
bursting  of  bubbles  at  the  sea  surface  ( Blanchard  and  Wbodcock,  1957).  The 
bubbles  are  produced  by  biological  activity,  chemical  reaction  and  breaking 
waves  (whitecaps)  .  Thus,  the  sea  surface  is  a  continuous  source  of  sea-salt 
aerosols  in  the  marine  boundary  layer.  At  wind  speeds  greater  than  about  3 
meters  per  second ,  whitecaps  are  the  primary  contributor  to  the  bursting  bub¬ 
bles.  These  surface  produced  aerosols  can  be  characterized  by  a  surface  flux 
spectrum,  Fg(r),  vhich  represents  the  volume  of  aerosol  per  particle  radius 
interval  produced  per  square  centimeter  of  ocean  surface  each  second  as  a  func¬ 
tion  of  aerosol  particle  radius.  This  quantity  is  a  function  of  wind  speed. 

The  continuous  production  of  sea-salt  aerosols  is  balanced  by  several 
removal  mechanisms .  Cne  obvious  mechanism  is  the  loss  of  particles  as  they 
fall  back  to  the  surface .  This  settling  under  gravity  is  called  "Stokes  fall¬ 
out"  and  is  characterized  by  the  Stokes  velocity.  The  particles  are  trans¬ 
ported  vertically  by  turbulence  in  the  marine  boundary  layer  and  maintained  at 
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a  nearly  ini  form  mixing  ratio  throughout  the  mixed  layer.  The  growth  of  the 
height  of  the  layer  constitutes  another  loss  mechanism  called  "entraiment" 
(Deardorff,  1976).  The  final  loss  mechanist  is  "rainout”  which  occurs  when 
the  particles  become  condensation  nuclei  in  the  formation  of  clouds. 

Given  the  surface  flux  spectrun  and  a  parameterization  of  the  removal 
mechanisms,  one  could  predict  evolutions  of  the  aerosol  density  spectrum 
(Ehirall  et  al,  1982).  Uh  fortunately,  the  surface  flux  spectrun  is  not  known. 
However,  we  can  reverse  the  process — that  is,  use  parameterize tions  of  the 
removal  processes  and  evolutions  of  the  aerosol  spectrun — -to  obtain  estimates 
of  the  surface  flux  spectrun.  This  paper  describes  the  calculation  of  Fg(r) 
from  data  taken  by  the  NFS  Environmental  physics  Group  during  the  Cooperative 
Experiment  for  West  Ooast  Oceanography  and  Meteorology  (CEWCCM-78)  and  the 
Joint  Air-Sea  Interaction  Experiment  ( JASIN) .  Details  about  the  experiments, 
equipment,  measurements,  and  analysis  are  given  by  Ehirall  (1981). 


2.  AEROSOL  SPECTRUM 


The  first  step  to  analyzing  the  evolution  of  the  aerosol  spectrin  is  to 
remove  those  variations  due  to  changes  in  the  ambient  relative  tumidity.  This 
is  done  by  transforming  each  spectrum  to  a  reference  relative  humidity  (RH  = 
00%  or  a  saturation  ratio  S  =  0.8) .  We  make  the  usual  assumption  that  the 
particle  radius  at  saturation  S  is  given  by  (Fitzgerald,  1975) 

rg  =  r-g(S)  (1) 

where  r  is  the  radius  at  S  =  0.8  and  g(S)  =  0.81  exp  [0. 066S/ (1. 058-S) ].  If 
we  have  a  measured  aerosol  number  density  spectrum,  n'(rs),  then  the  trans¬ 
formed  spectrum  at  standard  saturation,  n(r),  is: 

dn/dr  =  n(r)  =  n' (r-g(S) ).g(S)  (2) 

Since  we  prefer  to  work  with  the  volume  spectrum  V(r)  =  4/3wr  n(  r) ,  the  cor¬ 
responding  relationship  is: 

dV/dr  =  v(r)  =  V' (r-g(S) )/g2(S)  (3) 

Vfe  now  must  separate  the  measured  spectrum  into  its  two  basic  components : 

1)  The  background  continental  aerosol  of  nonlocal  origin,  V  ,  which  is 
present  above  and  within  the  marine  layer;  and 

2)  The  sea-salt  aerosol  locally  generated  at  the  sea  surface,  Vg,  which 
is  well  mixed  throughout  the  marine  layer  (h  is  the  height  of  the 
mixed  lajer)  .  Thus,  we  assume: 

V(r)  =  Vc(r)  +  Vs(r)  Z  <  h  (4a) 

V(r)  =  Vc(r)  Z  >  h  (4b) 

The  continental  aerosol  component  is  represented  with  a  Junge  type  dis¬ 
tribution  , 

Vc(r)  =  A/r  (5) 

where  A  is  the  continental  coefficient.  Since  there  is  very  little  sea-salt 
volume  production  for  r  <  0.3  ym,  we  can  use  the  snail  size  range  of  the 
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spectrum  to  calculate  the  continental  coefficient  (A  =  rV(r))  .  This  is  docu¬ 
mented  in  Fig .  1  where  ensemble  averages  of  volune  spectra  indicated  that  V(r) 
is  essentially  independent  of  wind  speed  at  r  =  0.1  pm.  At  low  wind  speeds, 
the  r  =  0. 1  pm  and  r  =  0.  3  wn  ensemble  average  spectra  are  described  quite 
accurately  by  Eqn  (5)  •  The  claim  that  A  is  an  accurate  index  of  continental 
influence  is  nicely  validated  by  comparison  with  atmospheric  ted on  activity 
(lars^n  et  al,  1979)  from  CEWCCM-78  (Fig.  2).  Thus,  the  locally  generated 
sea-salt  component  can  be  calculated  using: 

Vs(r)  =  V(r)  -  A/r  (6) 

The  rauainder  of  this  report  will  deal  exclusively  with  the  volune  spectra 
transformed  to  standard  humidity,  Vg(r)  .  We  write  this  as  Vs  to  simplify  the 
notation — with  the  understanding  that  Vg  refers  to  the  spectrum  at  RH  =  80%, 
with  radius,  r,  an  implied  variable. 

The  evolution  of  the  aerosol  spectrum  in  the  well  mixed  marine  boundary  is 
described  by  the  following  equation  (Fhirall  et  al,  1982), 

hdVs/dt  =  Fs  -  (We  +  Ws)  Vs  (7) 

( which  is  presented  here  without  further  discussion) 

where  Wg  is  the  entrainment  velocity  and  Wg  is  the  Stoles  velocity  (WU,  1979) 

Ws  =  1.57  x  10"2  (1  +  1.2(0.8l/g(S))3)r2  g2(S)  (8a) 

We  =  dh/dt  -  W  (8b) 

where  r  is  in  ym,  Wg  is  in  cm/ sec  and  W  is  the  mean  vertical  air  motion  due  to 

synoptic  scale  weather  (subsidence)  .  In  Bqn  (7 )  we  lave  left  out  the  cloud 

fomation  removal  mechanism  because  it  acts  on  a  much  longer  time  scale,  cne 

assumption  implicit  in  Eqn  (7)  is  that  the  convective  mixing  velocity,  W* 

(teimal  et  al,  1976),  is  much  larger  than  W  .  Since  W*  is  typically  1  n\/sec 

'  s 

while  W  is  cn  the  order  of  0. 1  m/sec  for  the  largest  particles  ( r  =  15  urn) 
s 


considered  here,  the  assumption  is  reasonable. 


dV/dr  (#fjm2/cm 
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3-  EUJX  SPECTRUM  ANALYSIS 

The  CEWCCM-78  data  have  been  used  to  evaluate  the  surface  flux  term  of  Bqn 
(7)  .  Since  the  other  terms  of  the  expression  were  measured,  one  simply  calcu¬ 
lates  F  using  a  rearrangement  of  the  terms.  Thus  the  final  expression  for 
s 

the  flux  calculation  is, 

Fs  *  MVs/dt  +  (We  +  WS)VS  (9) 

A  period  was  chosen  from  CEWOCM-78  where  all  data  were  available,  the 
synoptic  conditions  were  fairly  stable,  the  wind  speed  was  reasonably  constant 
and  a  good  mixed  layer  was  present  (Thble  I)  .  There  was  a  slowly  decreasing 
continental  influence  indicated  by  steady  NW  winds  (A  in  Table  I  and  Fig.  2). 

Application  of  Bqn  (9)  requires  a  value  for  the  entrairment  velocity.  If 
the  subsidence  velocity  is  known,  then  Wg  can  be  calculated  from  dh/dt  using 
Dgn  (8b)  .  Since  it  is  difficult  to  obtain  Vf  to  0.1  cny/s  accuracy,  it  was 
decided  to  estimate  Wg  empirically  fran  the  evolution  of  the  mixed-layer  tem¬ 
perature  and  water  vapor  density.  This  was  accomplished  by  simulating  the 
analysis  period  with  a  dynamic  mixed- layer  model.  After  initializing  the 
model  with  the  atmospheric  parameters  at  the  beginning  of  the  period,  the  en¬ 
trainment  velocity  was  adjusted  to  yield  best  agreement  with  the  temperature 
and  water  vapor  density  at  the  end  of  the  period.  The  mean  Wg  for  the  period 
was  0.35  an/s.  This  implies  W  =  -0.3  an/s  for  the  first  8  hours  and  W  =  0  for 
the  last  12  hours. 

The  scatter  and  uncertainty  of  the  aerosol  and  mixed  layer  depth  measure¬ 
ments  introduces  certain  trade-offs  between  averaging  times  and  statistical 
validity.  Ebr  short  averaging  times,  the  randan  scatter  of  the  individual 
terms  is  considerable ,  leading  to  large  variations  in  the  time  derivative 
terms.  By  increasing  the  averaging  time  to  4  hours,  the  F  " signal- to-noise" 
becomes  more  reasonable.  A  time  series  of  the  individual  terms  of  Efcjn  (9)  is 
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TABLE  I 


CEWCCM-78  meteorological  and  aerosol  data  for  the 
time  period  analyzed  for  this  report.  The  aerosol 
volume  spectra  density  dV/dr  is  given  at  each  radius 
r,  in  \m.  The  height  of  the  boundary  layer  was 
determined  fran  an  acoustic  sounder. 


DATE 

TIME 

U 

(nt/s) 

RH 

(%) 

h 

(m) 

A 

r=.  3 

r=.  8 

r=2  2 
( uni 

/an3) 

r=10 

r=15 

5/20 

1130 

7.4 

89.4 

420 

7.8 

15.1 

12.5 

9.79 

2.26 

0.801 

0.437 

5/20 

1230 

7.7 

87.6 

430 

5.2 

9.15 

12.5 

11.9 

2.08 

0.644 

0.327 

5/20 

1330 

8.2 

86.4 

430 

4.5 

8.39 

11.9 

12.2 

2.91 

1.11 

0.628 

5/20 

1430 

8.1 

82.1 

415 

4.7 

8.54 

10.6 

11.4 

2.24 

0.780 

0.420 

5/20 

1530 

8.3 

80.9 

437 

4.0 

7.96 

12.0 

11.2 

1.26 

0.296 

0.127 

5/20 

1630 

8.2 

81.1 

437 

4.0 

7.85 

11.9 

13.6 

2.47 

0.795 

0.411 

5/20 

1730 

9.0 

82.5 

443 

3.6 

7.26 

12.9 

13.4 

2.  55 

0.854 

0.455 

5/20 

1830 

10.0 

83.8 

430 

2.6 

6.54 

14.0 

13.2 

3.30 

1.37 

0-819 

5/20 

1930 

10.3 

84.3 

420 

2.4 

6.07 

13.8 

13.2 

3.71 

1.73 

1.11 

5/20 

2030 

9.8 

86.8 

360 

2.1 

6.54 

15.9 

14.2 

4.02 

1.63 

0.958 

5/20 

2130 

10.1 

87.7 

377 

1.9 

6.07 

14.9 

13.2 

3.34 

1.22 

0.678 

5/20 

2230 

9.5 

87.6 

397 

2.7 

6.82 

13.6 

12.8 

3.52 

1.42 

0.833 

5/20 

2330 

9.2 

87.6 

400 

2.7 

6.15 

12.9 

12.7 

3.46 

1.46 

0. 891 

5/21 

30 

9.6 

87.3 

420 

3.0 

6.14 

13.1 

11.8 

2.58 

0.993 

0.569 

5/21 

130 

9.2 

88.3 

490 

2.0 

5.69 

13.5 

12.9 

3.33 

1.26 

0.710 

5/21 

230 

8.3 

87.4 

520 

1.7 

5.29 

13.3 

13.4 

3.44 

1.39 

0.818 

5/21 

330 

8.2 

86.1 

560 

2.7 

6.11 

11.8 

10.9 

2-18 

0.  731 

0.390 

5/21 

430 

8.5 

83.1 

597 

2.2 

5.55 

11.6 

9.71 

1.78 

0.659 

0.369 

5/21 

530 

8.2 

81.3 

620 

2.0 

5.33 

10.9 

9. 75 

1.85 

0.732 

0-431 

5/21 

630 

7.9 

79.7 

680 

1.8 

5.01 

11.1 

10.7 

1.96 

0.678 

0.365 
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shown  in  Fig.  3  for  two  different  particle  sizes.  The  average  surface  pro¬ 
duction  flux  spectrun  is  shovn  in  Fig.  4.  This  flux  spectrun  applies  to  the 
average  surface  conditions  for  the  entire  20  hour  period  (average  wind  speed 
vas  about  9  n\/sec)  .  The  error  bars  represent  the  incertainty  in  the  mean 
estimate . 

The  relative  contributions  of  entrainment,  gravitational  fallout  and 

surface  production  are  nioely  illustrated  by  defining  an  equivalent  surface 

production  vertical  velocity,  W„  =  F_/V  so  that  in  equilibrium  (dV  /dt  =  0) 

r  s  s  S 

wF  =  We  +  Ws  (10) 

Ebr  the  average  conditions  fornd  in  the  analysis  period,  the  entrainment  and 
Stokes  terms  were  roughly  equal  at  a  particle  radius  of  5  um  (Fig.  5)  . 


Flux  (nm2/cm2/sec) 
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Figure  3.  Time  series  of  the  terms  of  Eqn  (9)  for  2  pm  radius  and  15  urn  radius 

particles.  The  dashed  line  is  the  average  contribution  of  the  surface 
production  term  and  should  be  the  sum  of  the  other  three  terms. 


Volume  Flux 
(|jm2/cm2/sec) 


Velocity  (cm /sec) 


Figure 


.  Ensemble  average  contribution  of  entrainment  (Wg) ,  Stokes 
fallout  (W  )  and  surface  flux  Wp  =  Fg/Vg  assuming  a  state 
of  dynamic  equilibrium. 
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4.  FLUX  SPECTRUM  AND  WIND  SPEED 


Because  of  the  reasonably  constant  wind  speeds  during  the  20  hour  CEWOCM- 
78  period,  we  were  able  to  improve  the  statistical  certainty  for  the  Fg  calcu¬ 
lation  by  combining  all  the  data  fran  the  period,  assuming  that  the  flux  would 
be  a  reasonable  representation  for  U  =  9  m/ sec  wind  speed.  Thus,  vre  now  have 
the  surface  flux  spectrum  at  a  single  wind  speed.  In  order  to  estimate  the 
flux  at  other  wind  speeds,  we  note  that  the  right  hand  side  of  Bqn  (10)  is 
nearly  independent  of  wind  speed  for  equilibrium  conditions.  Therefore,  the 
flux  at  one  wind  speed  can  be  related  to  the  flux  at  other  wind  speeds  if  the 
equilibrium  volume  spectra  are  toiovn : 

FS(UX)  =  fs(u2)  vs(u1)/vs(u2)  (11) 

Vfe  have  available  from  JASIN  a  large  set  of  ensemble  averages  of  aerosol 
volume  spectra  at  different  wind  speeds  (Fig.  6  and  Table  II).  It  is  a  simple 
matter  to  apply  this  data  to  Fiji  (11)  using  the  CEWCCM-78  aerosol  flux  and 
equilibrium  spectrum  to  generate  the  surface  volume  flux  spectra  as  a  function 
of  wind  speed  (Fig.  7  and  Table  III) .  At  a  wind  speed  of  6.0  m/sec,  the  flux 
can  be  integrated  in  radius  space  to  yield  a  total  sea-salt  dry  mass  produc- 
ticn  rate  of  4.8  mg/m  /day  which  compares  well  with  the  4.4  mg/m  /day  esti- 
mated  by  Hlanchard  (1963)  and  4.2  mg/m  /day  estimated  by  Kritz  and  Rancher 
(1980) . 


2  3 

TABLE  IX.  Bjuilibriun  sea-salt  aerosol  spectra  in  urn  /an 

as  a  function  of  radius  (mn)  and  wind  speed  (m/sec). 


r  ♦ 

0.8 

2 

5 

10 

15 

u 

6 

2-8 

3.4 

3 

0.3 

0.05 

9 

10 

10 

5 

1 

0.35 

11 

18 

25 

13 

6.4 

4.2 

13 

20 

30 

20 

15 

15 

15 

24 

33 

22 

22 

22 

18 

28 

35 

28 

28 

28 

* 

9 

9 

11 

2 

0.8 

0.4 

2  2 

TABLE  III.  Surface  sea-salt  aerosol  flux,  Fg(r),  (pm  /an  /sec) 
as  a  function  of  particle  radius  (pm)  and  wind 
speed  (n\/sec) . 


r  * 

\  u 

0.8 

2 

5 

10 

15 

6 

1.3 

1.1 

2.5 

1.0 

0.33 

9 

4.5 

3.1 

4.2 

3.3 

2.3 

11 

8.2 

7.7 

11 

21 

27 

13 

9.1 

9.2 

17 

49 

48 

15 

11 

10 

19 

72 

140 

18 

17 

11 

24 

92 

180 

* 

9 

4.1 

3.4 

1.7 

2.3 

2.6 

*  CEWCCM-78  Analysis  Period 
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5.  ^-EQUILIBRIUM 


Ebr  time  periods  of  a  few  hours,  the  aerosol  spectrun  may  not  be  in  a 
state  of  dynamic  equilibrium  If  we  rewrite  Bgn  (7)  in  the  following  form: 

dVs/dt  *  V,p  -  Fs/h 
where  the  time  constant,  t  ,  is 

Tp  =  V(We+Ws)  (13) 

then  we  see  the  analogy  of  aerosols  and  a  capacitor  charged  by  an  applied 
"voltage",  F_,  through  a  "resistance",  (W_  +  W  )_1.  In  this  analogy,  the 

S  6  S 

"capacitance"  is  h. 

The  response  time  of  the  aerosol  density  is  a  strong  fuicticn  of  particle 
2 

radius  because  W  -»  r  .  Vblues  of  t  for  S  =  0.8  and  h  =  400  m  are  given  in 
s  p  3 

Table  IV. 

TABLE  IV.  The  equilibrium  time  constant  for  aerosols 
at  S  =  0.8,  h  =  400  and  Wg  =  0.4  an/ sec. 

r,ym  0.5  1  5  10  15 

t  .  hours  28  22  11  3.5  0.5 

The  boundary  layer  mixing  time,  which  represents  the  time  required  for 
changes  in  aerosol  density  to  be  evenly  distributed  throughout  the  marine 
layer  is: 

Tm  =  h/W*  (14) 

Ebr  the  CEWOCM-78  analysis  period  (h  =  400  m  and  W*  =  0.6  m/sec)  we  find  rm  = 

0.16  hours.  Therefore,  short  term  variations  in  mixing  volume  (dh/dt)  will 

lead  to  changes  in  the  aerosol  density  because  the  production  response  time  is 

much  slower  than  the  mixing  time.  Thus,  for  time  periods  on  the  order  of  one 

hour,  changes  in  the  aerosol  density  (hdV_/dt)  will  be  highly  correlated  with 

the  mixing  volume  term  (V  dh/dt).  This  effect,  which  is  particularly  notice- 

s 

able  for  smaller  particles  (see  Fig.  3  and  Table  III),  is  shown  in  Elg.  8. 
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(oas/uio)  ^  tj 


-1.5  -1.0  -0.5  0  0.5  1.0  1.5 

-dh/dt  (cm/sec) 


Figure  8.  Changes  in  aerosol  spectral  density  (dV  /dt)  versus  changes 

in  mixing  volume  (-dh/dt).  This  graph  Illustrates  the  dominance 
of  these  terms  in  Eqn  (9)  for  short  time  periods. 
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